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. Abstract. We review our recent work on the classification of finite time singularities that 

£\J ' arise in isotropic universes. This scheme is based on the exploitation of the Bel Robinson 

energy in a cosmological setting. We comment on the relation between geodesic completeness 
and the Bel Robinson energy and present evidence that relates the divergence of the latter to 
^ ' the existence of closed trapped surfaces. 
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' 1. Introduction 

In we provided necessary conditions for the appearance of finite time singularities in isotropic 
universes based entirely on the behaviour of the Hubble parameter. These conditions provided 
us with a new classification of such singularities, a scheme which was further refined and 
expanded in j^j to include all known types of singularities and predict new ones. This second 
approach exploited the Bel Robinson energy, a kind of energy of the gravitational field projected 
in a sense to a slice in spacetime. Combining the quantities that constitute the Bel Robinson 
energy with the pressure and density of the matter fields through the field equations, we were 
thus led to a clear picture of how the matter fields influence the nature of the singularities. A 
complete classification was then obtained with the use of the Hubble parameter and the Bel 
Robinson energy. 

In the next Section, we introduce the Bel- Robinson energy and show that a closed Robertson- 
Walker (in brief RW) universe is geodesically complete if this energy is bounded. We present 
a summary of our classification scheme in Section 3, and in Section 4 we give necessary and 
sufficient conditions for the character of the resulting singularities. In the last Section, we 
discuss the existence of closed trapped surfaces in terms of the Bel-Robinson energy. 

2. Bel Robinson energy and completeness 

Consider a sliced spacetime with metric 

g = -N 2 {9 ) 2 + 9ij 9 l 9 j , 9° = dt, 6 i = dx l + ftdt, (1) 



where N = N(t, x l ) is the lapse function and (3 l (t, x 3 ) is the shift function, and the 2-covariant 
spatial electric and magnetic tensors 

= RiQj , 

D — —n n jjhklm 
J-^ij — ^'lihk'ljlm 1 ^ t 

Hij = 2^ lr lihkRoj' 

Bji = 2^~ lr ^ ihk ^j : 

where 77^ is the volume element of the space metric g. These four time-dependent space 
tensors comprise what is called a Bianchi field, (E, H, D,B), a very important frame field used 
to prove global in time results, cf. [3]. 

The Bel-Robinson energy at time t is given as the space integral 

B(t) = \l {\E\ 2 + \D\ 2 + \B\ 2 + \H\ 2 ) dp-n, (2) 

where by |X| 2 = g kl X^Xji we denote the spatial norm of the 2-covariant tensor X. In the 
following, we focus on a RW geometry filled with various forms of matter and with a metric 
given by 

ds 2 = -dt 2 + a 2 {t)da 2 , (3) 

where da 2 denotes the usual time-independent metric on the 3-slices of constant curvature k. 
For this spacetime, the norms of the magnetic parts \H\, \B\ are identically zero while \E\ and 
\D\, the norms of the electric parts, reduce to 



\E 



2 



3 (a/a) 2 and \D\ 2 = 3 ((a/a) 2 + k/a 2 ^ . (4) 



Therefore the Bel-Robinson energy becomes 

C 
~2 



B(t) = ^(\E\ 2 + \D\ 2 ), (5) 



where C is the constant. 

It is not difficult to show that in a closed RW universe such that \D\ is bounded above, H 
must be bounded above and the scale factor bounded below. Therefore H must be integrable 
and the spatial metric bounded below, that is such a universe is regularly hyperbolic. This 
in turn means that all hypotheses of the completeness theorem proved in |J| are satisfied and 
therefore such a universe is g— complete. It is also straightforward to see that the null energy 
condition is equivalent to the inequality \E\ < \D\, hence completeness is then accompanied 
with \E\ being bounded above. (For a flat spatial metric we would need to impose the regular 
hyperbolicity hypothesis in order to conclude completeness since the latter is independent from 
the boundedness of \D\.) 



3. Classification of singularities 

We are interested in finding the possible types of singularities that can arise in the evolution 
of a RW geometry. These types result from the possible behaviours exhibited by the different 
triplets consisting of the scale factor a, the Hubble expansion rate H and the Bel Robinson 
energy B. Assuming that the model has a finite time singularity at t = t s , the possible 
behaviours of the functions in the triplet (H, a, (\E\, \D\)) are classified as follows: 



51 H non-integrable on [ti,t] for every t > ti 

52 H — > 00 at t s > t\ 

53 H otherwise pathological 
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< CXD. 



The nature of a prescribed singularity is thus characterized completely by specifying the 
components in triplets of the form 

with the indices I taking their respective values as above. 

For fluid-filled RW models the various behaviours of the Bel-Robinson energy density can 
be written equivalently in terms of the density and pressure of the cosmological fluid: 

Bi ■v^ [i — ► cxd and \\i + 3p| cxd 

B 2 44> \i — > 00 and |/i + 3p| < cxd 

B3 44> jj, < cxd and |/i + 3p| — > 00 44> < cxd and \p\ — > cxd 

-B4 44> /i < 00 and |/i + 3p| < 00 44> < 00 and |p| < cxd. 

Following .2., we may gain further insight into the nature of any of the singularities 
expounded above by studying the asymptotic behaviour of the three functions that determine 
the character of the singularity. This is done by finding the relative strength, i.e., the limit as 
t — > t s of the relative rate of decay, for each pair of these functions. As an example compare 
the standard big bang singularities in the dust and radiation models in general relativity; even 
though both of singularities are of the type (Si,Ni,B\), they are really different since the 
radiation one has relative strength a « H « (\E\ <-> \D\) (where '<<', '~' mean that 
the relative limit is 0, 1 ore / 0, 1, with c a constant, respectively), whereas that in the dust 
models is characterized by the asymptotic behaviour a « H « (\E\ ~ 

4. Necessary and sufficient conditions for an (Si, Nj, B{) type of singularity 

Flat RW universes are clearly the simplest to analyze but in the long run one is interested 
to see whether the behaviours known to hold in the flat case continue to remain valid in 
more generic curved models. In this Section, we give necessary and sufficient conditions for 
the singularities in given flat RW models to continue to hold in more general universes in 
the sense of either having nonzero values of k or described by solutions in which some or all 
of the arbitrary constants remain arbitrary (particular or general solutions, cf. for this 
terminology). Usually the curvature term in a typical Friedman-type equation turns out to be 
subdominant compared to the density term or, in any case, does not alter the H behaviour. 
Some examples of what can happen in general relativistic models are shown in the following 
theorems. 



Theorem 1. Necessary and sufficient conditions for an (S2, N%, B±) singularity occurring at 
the finite future time t s in an isotropic universe filled with a fluid with equation of state p = wfi, 
are that w < — 1 and \p\ -^00 at t s . 

Theorem 2. A necessary and sufficient condition for an (£2, N2, B\) singularity at t s in an 
isotropic universe filled with a fluid with equation of state p + ^x = —B^, (3 > 1, is that /i — > 00 
at t s . 

Theorem 3. A necessary and sufficient condition for an (S3, N2, -B3) singularity at t s in an 
isotropic universe filled with a fluid with equation of state p + [i = — C(^o — z") -7 ! 7 > 0, is 
that fj, — > /io at t s . 

Theorem 4. A necessary and sufficient condition for an (Si, N%, B\) singularity at t\ in an 
open or flat isotropic universe filled with a fluid with equation of state p + fi = 7/x A ; 7 > and 
A < 1, is that fj, — > 00 at t\. 

Theorem 5. A necessary and sufficient condition for an (S\,N\,B\) singularity at t\ in an 
isotropic universe with a massless scalar field is that (ft — > 00 at t\. 

We can easily find examples of the behaviours described in the above theorems. The 
phantom cosmologies of [H] are described by Theorem the various solutions of jjj can 
be all accommodated by theorems (J2J and (j3J) respectively. The graduated inflationary models 
constructed by Barrow in |Hj when A = 3/4 belong to the situation of Theorem Q whereas the 
singularities of massless scalar field model of Ref. jH] are described by Theorem We shall 
meet the latter two models again in the next Section. The proofs of theorems 0-© and the 
details of the dark energy examples can be found in [2], our basic reference for all these results. 

One of the virtues of our classification scheme is that it can pick up even very mild 
singularities. The mildest in a sense type of singularity known to the authors was discovered 
in [7]. This is a general relativistic, flat RW universe filled with a fluid satisfying the equation 
of state 

and for (3 = 1/5 it admits the exact solution a = a e r8/3 . Then H = (8/3)r 5 / 3 , H = (40/9)r 2 / 3 
and H = (80/27)r" 1 / 3 . Thus as r — ► 0, a, a, d, H, H all remain finite whereas H becomes 
divergent. In this universe the Bel Robinson energy at an initial time, B(0), is finite whereas 
its time derivative is 



B(r) 



2\H + 2HH) + 4 ( A + H 2 ) f-^ + HH) 



a' \a* J \ a 2 J 



(6) 



and thus B(t) — > 00 at r — > 0. This is like a B4 singularity in our scheme. Since the derivative 
of the Bel-Robinson energy diverges, we may interpret this singularity geometrically as one in 
the velocities of the Bianchi (frame) field. At t s the Bianchi field encounters a cusp and its 
velocity diverges there. 

The character of the singularities can be expressed in terms of the functions describing the 
matter fields. We shall end this Section with a result of this sort. Note that the necessary 
conditions for singularities given in and based on the behaviour of the Hubble parameter 
alone can be rephrased in terms of the electric parts of Bel Robinson energy. We have the 
following result concerning the nature of the spacetime singularities in such models. 



Theorem 6. Necessary conditions for null and timelike geodesically incomplete, globally and 
regularly hyperbolic FRW universes are that at a finite time: 



B\ \E\ — > oo and \D\ —* oo 44> fi — > 00 and + 3p| — > oo 

-B2 < 00 and \D\ — > cxd 44> /i — > cxd and |yu + 3p| < cxd 

.B3 I £7| — ► 00 and \D\ < cxd 44> \i < cxd and |/i + 3p\ — > cxd 44> /i < cxd and |p| — > 00 

-B4 I £71 < cxd and |Z)| < cxd 44> < 00 and |^ + 3p| < cxd 44> < 00 and |p| < cxd. 



5. Incompleteness and Bel-Robinson energy 

It was recently shown by Ellis in JJ| that a RW space with scale factor a(t) admits a past 
closed trapped surface if the following condition is satisfied: 



a(t) > 



f(r) 



f(r) 



(7) 



with f(r) = sinr, r, sinhr for k = 1,0,-1 respectively. Recall that a closed trapped surface 
is a 2-surface with spherical topology such that both families of incoming and outgoing null 
geodesies orthogonal to the surface converge. Since our function \D\ can be written in the form 



\D\ 



V3_ 



d(t) 



f'(r) 
f(r) 



a(t) + 



f'(r) 

f{r)J ' P(r) 



+ 



(8) 



we see that the condition for the existence of a closed trapped surface becomes equivalent to 
the following inequality: 

We thus conclude that collapse singularities (as predicted by the existence of a trapped surface) 
are characterized by a divergent Bel Robinson energy. 

Consider for example the flat graduated inflationary model of [S] . This universe satisfies the 
null energy condition + p = 7/x 3 / 4 > 0, 7 > and the scale factor is given by 



a(t) = exp 



16 



3 3 / 2 7 2 t 



so that the inequality 



a(t) 



16 1 



3 3 / 2 7 2 t 2 6XP V 3 3 / 2 7 2 £, 

is satisfied for large r. We conclude that there exist past closed trapped surfaces in this model 
and by the singularity theorems we find that this model must be null geodesically incomplete. 
The nature of this singular behaviour is described by saying that it is of type (Si, N\,Bi). 

Another example is the flat RW filled with a massless scalar field given in [S] . This satisfies 
the null energy condition and has a past closed trapped surface: for large r we have 



16 



1 

> - 

r 



(10) 



(11) 



a(t) 



1 1 

3^73 > ~r 



(12) 



It is therefore null geodesically incomplete of type (5*i, Ni,B\). 

As a last example, consider the sudden singularity first introduced in jlOj which can arise 
in an RW model with scale factor given by 



<t) = ( j- 



1) + 1 



1 

ts 



(13) 



where 1 < n < 2 and < q < 1. We see that at t s , a(t) — > a s , d(t) — » a s and d(t) — > — oo. 
When fi >0 and p > 0, we find that there exists a future closed trapped surface (the analogous 
inequality of © is a(t) < -\ f \r) / f (r)|) since we have that the condition 

a(t) = l^'K " 1) + " i) n_1 < -I cot r| (14) 

is satisfied for r = 7r/2, a s < 1 and < t < t s . Therefore this model is timelike and null 
geodesically incomplete in the future with \E\ — > oo but \D\ < oo (type (S3, ./V2, -B3)). 

Acknowledgments 

This work was supported by the joint Greek Ministry of Education and European Union 
research grants 'Pythagoras' and 'Heracleitus' and this support is gratefully acknowledged. 

References 

[1] S. Cotsakis, I. Klaoudatou, Future singularities of isotropic cosmologies, J. Geom. Phys. 55 (2005) 306-315. 

[2] S. Cotsakis, I. Klaoudatou, Cosmological singularities and Bel-Robinson energy, arXiv: gr-qc/0604029 

[3] Y. Choquet-Bruhat, Constraints and Evolution in Cosmology, in Cosmological Crossroads, S. Cotsakis and 

E. Papantonopoulos, Eds. (Springer, 2002). 
[4] Y. Choquet-Bruhat, S. Cotsakis, Global hyperbolicity and completeness, J. Geom. Phys. 43 (2002) 345-350 

(arXiv: gr-qc/0201057 ); Y. Choquet-Bruhat, S. Cotsakis, Completeness theorems in general relativity, in 

Recent Developments in Gravity, Proceedings of the 10th Hellenic Relativity Conference, K. D. Kokkotas 

and N. Stergioulas eds., (World Scientific 2003), pp. 145-149. 
[5] S. Cotsakis and J. D. Barrow, The dominant balance at cosmological singularities, preprint, arXiv: 

gr-qc/0608137 

[6] P. F. Gonzalez-Diaz, Axiom phantom energy, Phys. Rev. D69 (2004) 063522 (arXiv: hep-th/0401082 ). 

[7] S. Nojiri, S. D. Odintsov and S. Tsujikawa, Properties of singularities in (phantom) dark energy universe, 

Phys. Rev. D71 (2005) 063004 (arXiv: hep-th/0501025 1. 
[8] J. D. Barrow, Graduated inflationary universes, Phys. Lett. B 235 40-43. 
[9] S. Foster, Scalar field cosmological models with hard potential walls, arXiv: gr-qc/9806113 
[10] J. D. Barrow, Sudden future singularities, Class. Quant. Grav. 21 (2004) L79-L82 (arXiv: |gr-qc/0403084 ). 
[11] G. F. R. Ellis, Closed Trapped Surfaces in Cosmology, Gen.Rel.Grav. 35 (2003) 1309-1319. 



